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Hardware problems are the most detrimental issues to channel estimates in 
wireless communication systems. Because of the enormous number of 
antennas at the base station (BS) in cellular massive multiple-input multiple- 
output (MIMO) systems and because one radio frequency (RF) chain per 
antenna is required, hardware impairments in such systems will be quite 
severe. Many research publications have used a quality-cost tradeoff to 
adjust for RF unit hardware issues. In this study, we have taken a different 
approach by reducing the error floor caused by impairments in the predicted 
channels. Here are two steps to remedy the problem. In phase 1, a single 
active user channel in a single cell was calculated statistically rather than 
parametrically. In phase 2, a convex optimization approach was used to 
regularize the estimated channel in phase 1 to reduce error and provide a 
robust channel estimate. The results of our proposed procedure are measured 
by the normalized minimum mean squared error (NMSE) versus a range 


from the effective signal-to-noise ratio, and it shows a significant reduction 
(nearly one order of magnitude) in the error floor as compared with the 
conventional one, especially at high signal-to-noise ratio (SNR) in the range 
of (20 dB-30 dB). Simulation results were extracted in MATLAB R2020a. 
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1. INTRODUCTION 

Massive multiple-input multiple-output (MIMO) systems are appreciated as an emerging key 
technology for future generations of wireless communication systems [1]—[8]. In most academic wireless 
communication research, an ideal hardware transceiver is assumed with a negligible hardware impact on the 
transmitted and received signals. In such a classical setup, both transmitter and receiver hardware 
components are characterized as a memoryless linear filter model. In other words, the transmitter correctly 
modulates its analog passband signal from the complex samples of its baseband signal, while the receiver 
demodulates and samples its received signal synchronously with the transmitter. 

However, the transceiver hardware in wireless communication systems is always non-ideal and the 
classical operations in ideal systems are not satisfied in practice [9]. In terms of power consumption and cost 
view, designing a massive MIMO system with a single radio frequency (RF) chain per antenna will lead to 
higher power consumption and expensive cost. Hence, a cost-quality trade-off policy should be followed, that 
is, the cost of realization of such systems will be M; times than that have been needed for the conventional 
systems with only one RF chain. Hence, lower-quality hardware components with a compensation procedure 
may be used for saving costs. 
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In effect, the aggregate non-linear impact of the various non-ideal hardware in the RF chain is called 
hardware impairments [10]-[15]. Because of the hardware impairments, there is a mismatch between the 
received signal and the baseband signal on the transceiver hardware ends. Different modeling systems with 
analog and/or digital algorithms such as in [16]-[18] and the references therein, have been issued to 
compensate for these nonlinearities problems and reduce their impact on the performance of the 
communication system. But none of these algorithms have completely fixed these problems in practice. In 
other words, residual impairments [19] still exist in many of these modeling systems raised either by 
inaccuracies in the design or by the destructive nature of some hardware impairments. 

In this paper, a single-user single-cell massive MIMO scenario has been considered under two 
impairment effects, the non-ideal hardware impairments, and the imperfection knowledge of the covariance 
matrices at the base station (BS) [20]—[22]. However, motivated by the fact that the channels estimated under 
such impairments have a large error floor, especially at high signal-to-noise ratio (SNR) [23], two phases of a 
procedure have been proposed here to estimate the channels under these conditions. In phase 1, the channels 
were statistically estimated based on many signal observations at the BS, exploiting the ergodicity of the 
channels and the low of the large number. In phase 2, to investigate a robust channel estimate with minimum 
error, a convex optimization method was applied to regularize the estimated channels in phase 1. Our results 
have been validated with conventional (parametric) channel estimation based on the local scattering model 
[24]. The results have shown a significant reduction in terms of floor error quantity at high SNR. 

The paper is organized as follows: in section 2, the concepts under the hood like residual 
impairments of the non-linear model, the uplink transmission model, and the channel estimation with 
hardware distortion have been studied in detail. Our proposed model has been explained in section 3, and the 
results and their analysis have been discussed in section 4. Finally, the paper is concluded in section 5. 


2. METHOD AND MODELS 
2.1. Residual impairments model (Bussgangs’ model) 

Bussgang theorem [25] suppose that if the complex Gaussian random variable like s ~ Cy (0, p) is 
fed into a distorting device with a non-linear memoryless function g(.), then the output of this device will no 
longer be Gaussian, that is if y = g(s), then y will be a non-Gaussian random variable, and both s and y can 
be cross-correlated as in (1). 


ya stn () 
where the term E{ys*}/p is a constant ratio that defines the correlation ratio between the output y and input 
variable s. The additive factor, 7, is a nonlinear term that defines the impact of the hardware device on the 
input s. Based on the Bussgang analysis, if we assumed E{|ys*|*}/E{s?} = p as a constant factor, and using 
this notation in (1), we can now generalize a consistent and analytical model to describe the aggregate effect 
of residual impairments on the RF signals. Hence, the input-output model that fully characterizes the 
detrimental impact of the residual hardware impairments in the communication system as in (2). 


y= /pstn (2) 


where s and y are assumed with equal power p. However, the power of the distortion term can be calculated 
from (2) as follows: {| |?} = E{lyl?}—p E{l sl} > {Inl} = (1 —p)p, that is n ~ Cy(0, (1 — p)p). 
Thus, the power of the distortion term is related to the input power p with a scaling factor (1 — p). That 
means, the additive distortion term is a power-dependent factor, in contrast to the classical additive receiver 
noise g? which is independent of the input power p. However, the constant factor, p € (0, 1], represents the 
hardware quality factor that can be used to measure the level of impairment of hardware devices at the 
transmitter and receiver as follows: when p = 1, the output will be the same as the input, that is, y = s, i.e., 
the distortion term ({| 7|?} = 0}), this is the case of ideal hardware. Otherwise, when p = 0, the distortion 
term will be 7 ~ Cy (0, p) which is the same Gaussian distribution as the input. The latter case represents the 
worst-case condition, i.e., the input s is completely distorted at the output of the non-linear device. 


2.2. Uplink transmission model 

In this section, we will study the aggregate impact of the hardware impairments on the RF signal 
from both the BS and terminal user (TU) transceivers, in addition to the propagation channel impact. In a 
multi-user scenario, assume the transmitted signal by an arbitrary TU k located in cell # is sj ~ Cy (0, jx). 


Bulletin of Electr Eng & Inf, Vol. 12, No. 1, February 2023: 215-223 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 Oo 217 


Following the model in (2), the complex Gaussian signal will be distorted by the non-ideal hardware of the 
TU device and then transmitted over the channel as yp” Sig +1j_ instead of s;,. The scaling factor 
p; “denotes the quality of the TU device, and for tractable notation, it is assumed to be the same for all TUs 
in the network. we additive factor nik determines the transmitter hardware distortion, which can be 
distributed as Mie ~ C,(0, (1 - pt” Pir). In general, the received planner array signals that reach the M; 
antennas at the BS # from all T in the network are as in (3). 


¥, = Loar pen: 1 hj, Ae Pi” Sei +n) (3) 


where L denotes the number of cells in the network, and Kọ refers to the total number of TUs in each cell in 
the network. However, for a given set of channel realizations {h3,} within an arbitrary coherence block, the 
analog signal ý; represents a complex Gaussian signal that has zero mean and a conditional correlation 
matrix which can be given as in (4). 


Ey, 74 {ne jt = Dhar DEL pee hd, (hé,)” a) 


However, since E {| Ver” Sei tne } = ppi. Therefore, it can once again apply the model in (2) to 
determine the a of hardware impairments at the receiver end (the BS). Thus, the signal y,; is replaced 
by J pes ¥, + n7 S, where p25 € (0,1] denotes the quality of the hardware of the BS # which is also assumed 
to be the same for all BSs, and n75 e C™ż is the hardware distortion term of the BS #. For analytical 
convenience, it has been assumed that the distortion terms between the different RF chains that attached to 
the M; receive antennas are independent. In general, the marginal signal y; E C™s that is received from the 
UL transmission at BS 7 can be modeled as in (5). 


y; = pF (Shas E hô; i( Pi” Sei + nv) +p +n; Š 


where n; is the White Gaussian (WG) noise that has been added at the BS 7 . Finally, the signal y; represents 
a complex baseband signal that will be used in the following channel estimation model. 


2.3. Channel estimation model under hardware distortion 
In this section, an uplink channel estimation model for the signal in (5) will be derived under non- 
ideal transceiver hardware based on the received signal Y? E C™+*Tp , To this end, let sp; = ./ppi Qo; E CP , 


represents the transmitted pilot sequence used by a TU iin cell #. Hence, the uplink transmission model in 
(5) over T, transmission instances will be given as in (6). 


Y? = /pP 5 (Sha Z rA h ( Pei PE” Pti + (nf?)")) + GP + Nz (©) 


where NP e C™4*Tp is the additive noise on the M; receive antennas with i.i.d elements and are distributed 
as ~ Cy-(0,07,), the term GPs e C™4*Tp is the pies hardware distortion matrix, each column of which 
takes the same distribution as n75 in (5), while 3 € C™ represents the transmitter distortion that contains 
Tp transmission realizations from no in (5). However, to estimate the channel of a TU i in cell £ at the BS #, 
the BS first correlates signal Y? with the pilot sequence pp; of the desired TU and yields the processed signal 
Yei = Ypi For an instant, to estimate the channel of the kth TU in cell ż at the BS #, the received signal 
will be written as in (7). 


Vise = Y Px = [Pix us Cpe PST h, + Èe, DEPik\G k) Pri pE V pë BiT hô; + 
OE a Ae Se 
De ae RE 


Desired pilot interfering pilots 


L UNT BS * P, 
Èj-1 ye 1 Veh AGH V Pj + G3” Pix + N; Pjk (7) 
Transmitter distortion Receiver distortion Noise 
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where T, results from multiplying PLP ik when 9}; = Pik: The set P,;, denotes for all users using pilot 


p 
: : : srs 2 š 
sequences as the pilot sequence @j, of the desired TU k in cell 7. Since Ilo; | = Tp . The last three terms in 
(7) will be distributed as follows: N? Wx ~ Cy (0, ObLTplm; i; CHa ~ Cy (0, t,(1 - pr’ psx): and 
BS . . . . . p 
G3" Pjk ~ Cy (0 m; Tp Dj {n}) for given channel realizations. It is seen from (7) that the processed signal y jik 
is affected by an aggregate distortion from the transmitter and receiver hardware of all signaling 
transmissions in the entire network. Since these distortions are almost surely non-orthogonal to their pilot 
sequence, this implies that there are pilot contaminations between every TUs in the network. However, 
according to the Bussgang theorem, the processed signal Vii in (7) represents a non-Gaussian signal. 


Therefore, a suboptimal estimator will be used here instead of the exact one. In other words, a linear 
minimum mean squared error (LMMSE) estimator with a little performance loss will be used instead of the 
optimal MMSE estimator. In general, based on the processed signal Yei = Y? Pei the linear LMMSE 


estimator [24] hf, can be given as in (8). 


hs, = y Pei pr pF RG YI Y pi (8) 


where R; e C™4%M} denotes the covariance matrix of h; and Y, refers to the overall effects from all 


interfering TU, the additive receiver noise, and hardware distortions caused by pilot signaling of all TUs in 
the entire network, it is written as in (9). 


, ; e 
w= (Zene ew Pei PE” Pe Ri, + Stim, Pa Dida Pa- 

A K =f 
pr” )pPs Re, + Eiza Lica Pailt — PPS) Dye ) O) 


where 
Dp, = diag (Ril ay [Redlu, a (10) 


represents a diagonal matrix that contains diagonal elements from Rees However, the estimator that gives a 
minimum mean squared error minimum mean squared error (MSE) E{||hé, = i, | } is the best estimator that 
can be used, the difference hij; - hi, represents the estimation error he: that has a correlation matrix 


f A E RE, 
ci, = E {f4 (f4,) } given as in (11). 
Co, = Ro, — Pei PE” Pr Tp RG Pii Rg (11) 
while the correlation matrix of the linear estimator he calculated as in (12). 
{hé, (64) }= ti Mei = Pei Pt Pr Tp eit pipi (12) 


It is noticed from (8) that the correlation matrix of the linear estimator i is fully dependent on the 
second-order statistics of all TUs in the network, which means that the BS can determine the LMMSE 
estimator ie only when it knows the correlation matrices RŹ; and wi including the scaling quantity 
Pei Pi are However, the BS in practice has imperfect or no prior knowledge of this statistical 
information and therefore should be statistically estimated based on the random matrix theory. 


2.4. Hardware distortions impact 

In this section, a single cell single active TU scenario has been tested to show how hardware 
impairments impact the estimated channel in section 2.3, the channel between the TU and the BS is denoted 
by h ~ Cy-(Oy, R), where R E€ C”*™) is the semi-definite covariance matrix. Following (11), the correlation 
matrix of the error for a single user scenario is given as in (13). 


C =R- p pi’ pSt, RYR (13) 
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Where 


=1 
Y = (p (1 + p7” (tp — 1)) pPSR + pd — p25)Dp + oĝ,1m) 14) 


The Y expression in (14) shows the distortion effect of the transceiver’s hardware in the desired TU 
and the BS, which is different from the case when having ideal hardware i.e. p7” = p85 = 1. Sub. (14) in 
(13) and take the limits when p goes to infinity the correlation matrix of the error becomes as in (15). 


_p_ pftter (p21) p , (1-08) | \~ 
C=R R( pu, Rtas De) R (15) 


For a special case of R = Bly, the correlation matrix of the error in (15) simplifies to: 


_ _ g2 _ _Ba-p{" prs) 
C = Bim 1+ p14 (tp -1) (1-85) Iu = 1+p1 0 pPS (tp -1) (16) 
TU 

t Tp 


p pTU pBStp 

The last expression in front of Iy represents the aggregate error floor term of the transceiver’s 
hardware at the BS and the TU device. As seen in (16), the error floor is constrained by the hardware quality 
pi” and p25, the pilot length Tp , and the average gain of the channel B. Hence, one can use these parameters 
to optimize the error floor at the BS. 

The normalized mean squared error normalized minimum mean squared error 
(NMSE) = tr(C)/tr(R) in Figure 1 has been used in [24] to show how the hardware distortion affects the 
estimated channel when using “the local scattering model” with angular standard deviation (ASD) = 10° and 
a uniform distributed nominal angle. NMSE curves are averaged over a range from the effective SNR for 


Tp = 10 lengths and with different hardware constants, as illustrated in Figure 1. 


10° F 3 
10 F 
pt” = př° = 0.95 
2 < 
10° 
w 
Q 
= 
2 . 
10° F aa 
104 F ie 7 
Impairments: LMMSE estimator | 
oie Ideal hardware: MMSE estimator | 
10° 1 1 f i 1 J 
-10 -5 te) 5 10 15 20 25 30 


Effective SNR [dB] 


Figure 1. The normalized MSE of the conventional channel estimation using the local scattering model in 
[24] with equal hardware quality at the transceiver ends 


3. PROPOSED APPROACH 

In motivation with the fact that the BS can estimate the channel hź; only when it knows the 
correlation matrices RS, and wi which in practice are unknown at the BS. Therefore, in this section, the 
methods of two phases have been suggested to estimate the channels under such consideration. In phase 1, 
the channel hi, statistically estimated based on many observations of the signals in (6) exploiting the channel 
ergodicity characteristics and the low of the large number. A sample covariance matrix from the random 
matrix theory of the desired TU channel would be estimated first, then a convex optimization method will be 


applied to regularize the covariance matrices in phase 1. In the present work, two imperfection conditions, 
the hardware distortions, and the unknown statistics of the channel at the BS have been jointly considered. 
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3.1. Phase 1: sample covariance matrix estimation 
3.1.1. P% Sample) estimation 
In this part, the BS needs to see many observations from Viei as Yieiltl, Ge [Nw], where Ny 


refers to the total number of signal observations, assumed the coherence time is large enough to include many 
channel realizations, then the sample matrix can then be given as in (17): 


Di (sample) 


1 oN H 
Pi Np Un=1 Yei [n] (Pa [n]) (17) 
and according to the large number low, i.e., when Ny — ©, then as in (18). 


1 
Nw 


nVa (al) [A a8) 


which is a diagonal element from the sample covariance matrix in (17). 


3.1.2. R4, Sample) estimation 

In this part, the same approach above can be followed with two stages of a procedure as in [21]. In 
the first stage, the BS first estimate the wi (sample) that determined for all TUs in the network, including the 
desired TU. In the second stage, the BS gives extra pilots to the interfering users and used their observations 
to estimate RS, (sample) as in (19). 


aj; (sample) ~ ; (sample) ah (sample) 
A — j A 
Re; = P; — Yii -k (19) 
ay 7 (sample), ; ay (sample), i ; 
where W;, is the already computed matrix and P5; _, is the sample covariance matrix of the 


interfering TUs with additional pilots. 


3.2. Phase 2: convex optimization of sample covariance matrices 

Due to the estimation error in all elements of the sample covariance matrices estimated in phase 1; a 
convex optimization method will be used here to obtain a better estimate, where the statistical sample 
covariance matrices are regularized with its diagonal matrices to reduce the error as follows: 


ol ~ ; (sample) ~ ; (diagonal) 
Pc) = c Pi, +1- c)? (20) 
Pp. a; (sample) a; (diagonal) 
Ri; (a) = afb; + (1—-a)Ri, (21) 


where c and a denote the regularization factors. The diagonal elements of the Rs, (a) and pi (c) are the 
. sj (diagonal) Di (diagonal) . g 7 . , 
same as in Rj, and P3; respectively, while the elements in the off-diagonal locations are 


regularized corresponding to c and «æ € [0,1]. This might be important for the estimating model in (8) since it 
uses the covariance matrices to estimate the channel Hiss Thus, (8) can be rewritten in terms of the 
regularized covariance matrices as in (22). 


(a,c) = Vpr pr” PPR (a) Pf, (c) yha (22) 


Then, the NMSE = tr(C)/tr(R) for a single-TU single-cell scenario can be applied to compare the 
obtained results with the standard results in Figure 1. 


4. RESULT AND ANALYSIS 

For performance evaluation, the average NMSE curves in Figure 1 of the local scattering model have been 
considered as the standard results [24] that can be used to validate our results in Figures 2 and 3. It has been assumed 
that the NMSE curves in Figure 1 are averaged over different nominal angles of the TU in the cell (0° — 360°) with 
Gaussian distribution and angular standard deviation (ASD = 10°) using M = 100 antennas at the BS. Figure 1 
illustrates the NMSE curves as a function of the effective SNR (p/af,) with a pilot sequence Tp = 10 and different 
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hardware constants as follows, p7” = p85 = 0.99 and p7” = p85 = 0.95 in case of hardware distortions, and 
pi’ = pPS = 1 in case of ideal hardware. Figure 1 emphasizes that there is an error floor in the estimated channels, 
particularly in the range of (20 dB and 30 dB) SNR, while it is small at low SNR values. On the other hand, 
hardware impairments show a substantial impact on the estimated channels when the quality of hardware 
components decreases, particularly at p?” = p25 = 0.95. In the present work, an optimized channel estimation 
method was proposed to reduce the error floor in the NMSE curves at high SNR. We first take advantage of the 
random matrix theory to estimate the channels statistically based on the sample covariance matrix method under 
imperfect statistical information, then a convex optimization procedure has been applied to regularize the estimated 
sample covariance matrices of the estimated channels. Our results illustrated in Figure 2 show the average NMSE 
curves for the estimated channels under the joint impact of hardware impairments and imperfection knowledge of 
the channels at the BS. It shows how the practical channels are closely aligned to the standard NMSE curves of the 
local scattering model in Figure 1. On the other hand, Figure 3 shows the optimized NMSE curves when the convex 
optimization procedure was applied to the sample covariance matrices in Figure 2. As a result, the error floor in 
Figure 3 is significantly reduced when the hardware quality p7” = p2* = 0.99, while it is slightly better than 
before in the case of pł” = p®° = 0.95. On the other hand, it is noticed that the NMSE curve when 
pi’ = pēs = 0.99 is closely aligned to the ideal hardware case of p7” = p25 = 1, which means that less quality 
hardware components can be used now at the BS with only a little performance loss in the channel estimation. 


NMSE 


Impairments: LMMSE estimator bs 

al Impairments: Statistical LMMSE estimator 

= = = ‘Ideal hardware: MMSE estimator 

40° L 1 1 L Í "i 1 
-10 -5 0 5 10 15 20 25 30 


Effective SNR [dB] 


Figure 2. The normalized MSE of the statistical channel estimation using the sample covariance matrix in 
phase 1 with equal hardware quality at the transceiver ends 
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Figure 3. The normalized MSE of the optimized channel estimation using the convex optimization method in 
phase 2 with equal hardware quality at the transceiver ends 
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5. CONCLUSION 

In this research, a channel estimate approach for large MIMO systems was suggested under two 
situations of impairment: the non-ideal hardware impairments at the transceiver ends of the transmitter and 
receiver, and the impairments caused by the BS's insufficient statistical information. We estimated the 
channels jointly under these conditions, regularized them, then compared them to the typical result from the 
local scattering model. In terms of hardware quality constants, the results show that our proposed technique 
beats the usual one in terms of lowering the error floor (almost by an order of magnitude). 
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